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■■flrtnermokinetic SI engine simulation was used to study the effects of simple 
NO x control techniques on performance and emissions of a methanol fueled 
engine. As part of this simulation, a ring crevice storage model was formu- 
lated to predict UBF emissions. The study included spark retard, two methods 
of compression ratio increase and EGR. The study concludes that use of EGR 
in high turbulence, high compression engines will both maximize power and 
thermal efficiency while minimizing harmful exhaust pollutants. 
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Alcohol fuels are rapidly gaining popularity as excellent alternatives to 
petroleum fuels throughout the world. Studies in our laboratories have shown 
alcohols to be superior to gasoline in performance, emissions, health and 
safety [1,2].* Continuing research and vehicle fleet studies define problem 
areas which need to be solved before alcohol fuels are ready for public use. 

A united world effort will continue to reduce these problem areas such that 
alcohols can be safely introduced into the marketplace without concern for 
larger problems downstream. 



While alcohol fuels are generally clean burning, air pollution problems exist- 
ing in major cities throughout the world will require emissions controls on 
all combustion sources. Since oxides of nitrogen (NOx) emissions are funda- 
mental in the formation of photochemical smog and are the hardest emission to 
control, a review of simple NOx control techniques is presented in this paper. 
With this study, a simple method of N0 X control can be defined without large 
penalties in performance or exhaust emissions. 

The method used in this paper to study these N0>, control techniques is a fun- 
damental spark ignition (SI) engine computer model. This approach is used 
since the interaction between engine design and operating variables is 
sufficiently complex that a fundamental model is essential in interpreting 
performance and emissions results from a spark ignition engine. A fundamental 
model employs current Knowledge and semi -empirical formulations of combustion 
and pollutant formation processes occurring during an SI engine cycle to 
predict, in a reasonably expedient and low-cost fashion, the results of various 
changes in engine design and operating variables. The results from this kind 
of study can be used to guide experimentalists in interpreting their results 
and engine designers in screening various concepts in engine design prior to 
expensive hardware development and testing. 

THE COMPUTER MODEL 

The computer model used in this study has been under continuous development 
for the past ten years. During this time, various highlights of this model 
have been presented in the literature [3,4,5]. The model is a combination of 


*Numbers in brackets [] denotes references listed at end of paper. 
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two submodels, the kinetic Otto cycle submodel and the kinetic exhaust sub- 
model. The kinetic Otto cycle submodel includes both thermodynamic and chem- 
ical kinetic considerations of processes occurring during the intake, com- 
pression, combustion and expansion phases of the actual SI engine cycle. Its 
details are published in Reference 3. Modification of the model to incorpor- 
ate squish chamber cylinder head geometry is outlined in Reference 4. This 
model uses the detailed high-temperature chemical kinetic reaction mechanism 
listed in Table I to define methanol oxidation, and carbon monoxide (CO) and 
N0 X pollutant emissions formation which occur during the combustion and expan- 
sion phases of the engine cycle. 

The kinetic exhaust submodel is a one-dimensional fluid mechanics model of the 
exhaust port and manifold, and applies a detailed low- temperature kinetic 
reaction mechanism to predict unburned fuel (UBF) oxidation and aldehyde (ALD) 
formation during the exhaust event. The details of this model are discussed 
in Reference 5. The low- temperature reaction set uses all the reactions in 
the high temperature reaction set; however, reaction rates for reactions 24 
through 28 have been replaced with those listed in Table II to better predict 
aldehyde oxidation at low temperatures and two additional N0 X reactions have 
been added to better predict N0 2 formation during the exhaust stroke (also 
listed in Table II). 


Further refinements to the kinetic exhaust submodel result from recently pub- 
lished papers on hydrocarbon quench layer content [6,7,8]. These studies 
suggest that under normal operating conditions, wall quenching is not an 
important source of exhaust hydrocarbons. This is due to rapid diffusion and 
burn-up of quenched hydrocarbons on a millisecond time scale without formation 
of significant aldehydes [7]. These studies believe that exhaust hydrocarbons 
result from ring crevice storage and absorption-desorption of hydrocarbons in 
the lubricating oil. Earlier studies by Wentworth have shown up to 74% re- 
ductions in exhaust hydrocarbon concentrations by eliminating the ring crevice 
[9]. Thus, a submodel utilizing ring crevice storage was formulated to calcu- 
late input hydrocarbon concentrations for the kinetic exhaust model. This 
ring crevice storage model assumes that half of the unburned fuel existing in 
the ring crevice at the end of combustion exits the crevice during the expan- 
sion stroke and is subsequently burned up on the bulk cylinder gases before 
the exhaust valve opens. It is also assumed that 10% of the unburned fuel in 
the crevice does not exit the engine cylinder during the exhaust event. The 
remaining 40% exits the cylinder under the following assumptions. One-half 
of the exiting unburned fuel is evenly entrained in the bulk cylinder flow 
during the time 90% of the cylinder charge (minus exhaust residual) is ex- 
hausted from the cylinder. The remaining half exits the cylinder with the 
remaining 10% of the bulk cylinder charge but at five times the concentration 
level. These assumptions result in exiting hydrocarbon concentration profiles 
found in experimental studies [10,11]. 


The method above was used in calculating input hydrocarbon levels for the ex- 
haust model in all but the exhaust gas recirculation (EGR) study. In this 
part of the study, additional sources of unburned fuel emissions needed to be 
considered. One such source resulted from reducing manifold vacuum as EGR 
was added to maintain the same volumetric efficiency (nv)- This tends to 
increase fuel flow through the cylinder during the valve overlap period due 
to scavenging. Also as EGR is added, there is a general degradation from good 
burns to slow burns [12], thereby increasing unburned fuel emissions. Calcu- 
lations using only ring crevice volume resulted in decreasing input hydro- 
carbon concentrations with increasing EGR. Thus for this part of the study, 
the input hydrocarbon concentration was calculated for the no EGR case using 
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the ring crevice model and then used for all EGR addition cases as well. 

Run times for the entire engine cycle simulation were typically 1100 seconds 
on a CDC 7600 computer. This provided the best compromise between computer 
time and numerical accuracy. 

MODEL TUNING AND VERIFICATION 

In order to tune the computer model's semi -empirical equations, experimental 
data was taken from the four-cylinder 2.3 liter Ford Pinto engine used by the 
University of Santa Clara in methanol characterization studies under federally 
funded grants [1,2]. The specifications for this engine are given in Table III. 
Cylinder number one of this engine was instrumented with a model 601 B Kistler 
pressure transducer mounted in a model 640 spark plug adapter and connected 
to a model 504E charge amplifier and a Tektronix model 5011 storage oscilli- 
scope. Exhaust emissions were taken at the cylinder's exhaust port and 
analyzed. 

The corrected pressure trace obtained from the engine was then used to deter- 
mine coefficients for both the flame speed and heat transfer equations. As 
shown in Figure 1, good agreement was reached between the computer model and 
actual pressure diagram in combustion rate; however, the model tended to pre- 
dict higher peak pressures than the actual pressure trace. 

Two test points are compared against the model's predictions in Table IV. The 
first (<t> = 0.9, CR = 8.44)* represents the test conditions at which the cylin- 
der pressure trace shown in Figure 1 was taken. The second ($ = 1.0, CR = 9.0) 
represents the base engine operating condition used in the parametric study. 
Both comparisions show good agreement in performance and emissions predictions 
with the actual experimental data. Predicted CO emissions, however, were much 
lower than actual values. 

SPARK RETARD EFFECTS 

In a spark ignition engine, spark timing has a large influence on the combus- 
tion process and, thus, on power, thermal efficiency (nth) and emissions. As 
shown in Figure 2, N0 X emissions decrease rapidly as the spark is retarded 
from MBT (minimum spark advance for best torque) due to burning later in the 
cycle and, thus, at lower peak temperatures and pressures. CO emissions 
tended to decrease with spark retard due to higher temperatures during the 
expansion stroke. Unburned fuel emissions also tend to decrease with spark 
retard due to higher exhaust temperatures. Aldehydes, however, peak slightly 
retarded of MBT, then decrease with more retard. While the increasing exhaust 
temperatures tend to react more unburned fuel into aldehydes as the spark is 
retarded, the amount of unburned fuel in the ring crevice decreases with spark 
retard due to lower peak pressures. The latter becomes the overriding effect 
slighly retarded of MBT, thereby causing the aldehyde emissions to drop as the 
spark is retarded further. Pischinger [19] found similar results for lean 
mixtures, but found a steady increase in aldehyde emissions with spark retard 
for rich mixtures. 

COMPRESSION RATIO EFFECTS 

Due to the high octane rating of methanol, higher compression ratios can be 


* $ denotes fuel -air equivalence ratio while CR denotes compression ratio. 
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used without fear of engine knock. In conventional engines, squish chamber 
combustion chamber geometry is used to minimize the pressure rise rate during 
both the start and end of combustion. Reducing the pressure rise rate during 
the start of combustion, which occurs during the compression stroke, minimizes 
the negative work. Reducing the pressure rise rate at the end of combustion 
minimizes the possibility of knock by minimizing end gas heating. 

In this study, two methods of raising compression ratio are examined. The 
first method maintains the same open-to-squish chamber volume ratio at top 
center (TC) while compression ratio is raised. This maintains the same 
relative squish velocities with varying compression ratio and would have 
relatively the same effect in power and emissions as raising compression ratio 
in a disc chamber engine. This method is best accomplished by piston crown 
redesign. The second method maintains the same squish chamber clearance 
height while compression ratio is raised. Thus, compression ratio is raised 
by lowering only the open chamber height. This results In decreased squish 
velocities as compression ratio is raised (in comparison to method one), 
thereby reducing overall in-cylinder turbulence levels. This method is nor- 
mally accomplished by milling the cylinder head. 

With the first method, N0 X increases 29.6% by increasing the compression ratio 
from 9:1 to 15:1 due to higher peak pressures and temperatures as shown in 
Figure 3. Unburned fuel emissions rose 108.5% for the same comparison due to 
increased ring crevice storage and decreasing exhaust temperatures. Aldehyde 
emissions, however, decreased 43.9% due to lower exhaust temperatures , thereby 
reacting less unburned fuel during the displacement flow period of the exhaust 
stroke when aldehyde emissions are generally formed [5]. CO emissions rose 
31.5% for the same comparison while both power and thermal efficiency rose 
14.7%. 

In the second method, the reduced turbulence levels and burning rate tended to 
limit the increase in NOx and CO emissions as shown in Figure 4. By raising 
the compression ratio from 3:1 to 15:1 by this method, N0 X emissions rose 
21.8%, UBF emissions rose 109.8% and CO emissions rose 25.3%. Aldehyde 
emissions, however, decreased only 19.3%. The trends of aldehyde emissions 
with compression ratio shown in Figure 4 are similar to the trends found by 
Pischinger [19] using a similar method to raise compression ratio. Power and 
thermal efficiency rose 15% for the same compression ratio increase. 

EGR EFFECTS 

In this study volumetric efficiency was maintained as EGR was added. A throt- 
tled condition was chosen for this study to compensate for the increased 
cylinder charge due to EGR addition. Furthermore, the inlet charge tempera- 
ture was allowed to rise with EGR addition due to the mixing of the hot exhaust 
(600°K) with the cool incoming charge (289°K). This raised inlet mixture 
temperature substantially from 289°K for the no EGR case to 392°K for 30% EGR 
addition at the 65% volumetric efficiency condition. The decrease in inlet 
temperature as EGR was added more than compensated for the increased burning 
velocities, thereby raising both power and thermal efficiency as shown in 
Figure 5. Both power and thermal efficiency rose 8.5% with 30% EGR over the 
no EGR case. N0 X emissions dropped substantially showing a 97.2% reduction 
for the above comparison due to lower peak cylinder temperatures and pressures. 
These results are consistent with other investigators [20,21]. Unburned 
fuel emissions increased 55.8% for the above comparison due to decreasing 
exhaust temperatures and oxygen concentration. This matches trends found in 
our laboratory [21]. These factors, however, reduced aldehyde emissions 88.3%. 
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CO emissions dropped 70.7 % for the same comparison. The interesting point is 
that as the exhaust temperature decreases, more of the NO is converted to 
N0 2 . At the no EGR case, only 8.2% of the N0 X is N0 2 , while of the 30% EGR 
case, 49% of the N0 X is N0 2 . Thus, EGR Is an effective method of reducing 
N0 X emissions without large increases in UBF emissions while also substan- 
tially reducing CO and aldehyde emissions. 

CONCLUSIONS 

A thermokinetic computer model was used to study the effects of simple N0 X 
control techniques on performance and exhaust emissions. A ring crevice 
storage model was also formulated for predicting unburned fuel emissions. The 
study showed spark retard to be an effective method of reducing N0 X , CO and 
UBF emissions at the expense of power, thermal efficiency and aldehyde emis- 
sions. It was also shown that by raising compression ratio in a squish 
chamber engine, the rise of N0 X , and CO emissions could be limited, but not 
without an aldehyde emission penalty. Exhaust gas recirculation, however, 
was shown to be a very effective method for reducing N0 X , CO and aldehyde 
emissions with only a slight UBF emission penalty. It is the opinion of the 
authors, that use of EGR in high compression methanol fueled engines with 
increased turbulence to reduce cycle-to-cycle variatiors will both maximize 
power and thermal efficiency while minimizing harmful exhaust emissions. 
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Table III: Engine Specifications 


Engine Type 

Squish Chamber ONV 

Bore 

9.60 cm 

Stroke 

7.94 cm 

Connecting Rod 

13.22 cm 

Intake Valve Diam/tift 

4,41 an/1.016 cm 

Opens 

60° BTC 

Closes 

100° BTC 

Exhaust Valve Diam/Uft 

3.81 cm/1.016 cm 

Opens 

’00° ATC 

Closes 

60° ATC 

Fuel Type 

Methanol 

Intake Temperature 

284®K 

EGR Temperature 

600° K 

Ring Crevice Volume 

0.9072 cc 
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Table IV: 



* 4 * i.* a u** 

CNAN* ANGLE DEGREES (ATO 


Performance and Emissions Compa *isons at 
Match Points < WOT & 2000 RPN) 



ENGINE 

DATA 

MODEL 

ENGINE 

DATA 

mon 

Performance: 





♦ 

0.90 

0.90 

1.00 

1.00 

CR 

6.44 

8.44 

9.00 

9.00 

SA (-'BTC) 

34 

34 

?S 

25 

IMP (MO 

9.84 

9.76 

10.70 

10.8) 

n , h 

39. SO 

39 .ro 

39.20 

39.00 

Emissions : 





URf (PPM) 

4 f '9 

41 S 

S3S 

520 

ALOE mots (PPM) 

w; 

m 

148 

140 

CO (■) 

o.o; 

0.03 

o.rs 

0.09 

NO (PPfU 

?4U 

?3b6 

I9?S 

1894 

H0 % (PPM) 

?SS4 

7488 

2042 

I9S0 


I I CURE 1 1 COM! * AM I GON 
PRESSURE 


UP PQMPUT i R PREDICTED 
RAQ 10 ACTUAL DATA 


v * 






ie “T 

2 ^ » - 

sE 

<—v U. 

g Ui 30 L 
8.20 r 


HFTiri 

1.00 

9.88 CR j 

9.09 CR 

2008 RPR 

2M RPR 

N9T 

MBT 

%• 79.21 J 

V 79.2* 

METHANOL j 

i 1 

METHANOL 

t 12 3 fs *s r 


- 

p 

s li. 

a*s» is 



-ie *ai ie ra 

SPARK RETARD (DECREES) 



It 13 IS 

COMPRESSION RATIO 


FIGURE 2* SPARK RETARD EFFECTS 


^ - 1.00 
9.89 CR 
2000 RPR 
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FIGURE 3i COMPRESSION RATIO EFFECTS 
CONSTANT OPEN-TO-SQUISH VOLUME RATIO 


1.00 
9.88 CR 
2280 K . 
MBT 

V 651 



S& 35 



ALO 


O X 

UJ CL 
Z CL 

5 w 608 


9 11 13 lb 

COMPRESSION RATIO 

FIGURE 4i COMPRESSION RATIO EFFECTS 
CONSTANT SQUISH CLEARANCE HEIGHT 


n a 

§ ~ 588 



< 18 28 30 

ECR (X) 

FIGURE 5» EGR EFFECTS 
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